Magnetic Field Estimation at and beyond 1/N Scaling via an Effective Nonlinearity 
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We provide evidence, based on direct simulation of the quantum Fisher information, that 1/N 
scaling of the sensitivity with the number of atoms N in an atomic magnetometer can be surpassed 
by double-passing a far-detuned laser through the atomic system during Larmor precession. Fur- 
thermore, we predict that for N ~S> 1, the proposed double-pass atomic magnetometer can essentially 
achieve 1/N scaling without requiring any appreciable amount of entanglement. 
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PACS numbers: 07.55.Ge, 32.80.Pj, 33.55.Fi, 41.20.Gz 

Introduction — The strength of a magnetic field B is often 
determined by observingLarmor precession in a spin- 
polarized atomic vapor [l|, More generally, a sys- 
tem of N atoms can be prepared into a known initial 
state po and let to evolve to the final state p T {B) = 
U T (B)p Ul(B) at time t — t under dynamics Ut{B) pa- 
rameterized by the field B. After (or preferably during) 
the evolution, the atoms can be measured to infer the 
field strength from the dynamics using the methods of 
quantum parameter estimation theory [j, 0, 0, E|- 

For precise measurements, the uncertainty SB in the 
estimated value of the field strength B is dominated by 
quantum fluctuations in the measurements performed 
on the atoms. The quantum Cramer-Rao inequality 
0, 0] places an information-theoretic lower bound on the 
(units-corrected mean-square) uncertainty in terms of the 
quantum Fisher information It{B) 



SB T = 



B T 



\d(B T )/dB\ 




(1) 



which holds for any estimator used to determine B. For 
pure states, the quantum Fisher information is given 
by the expectation value 1t{B) = tr[£f(B)p t (B)] of the 
symmetric logarithmic derivative operator 



£ t (B) = 2dp t (B)/dB, 



(2) 



which characterizes the sensitivity of the time-evolved 
state Pt(B) to variations in the value of the parameter 
B. Subsequently, £t(B) can be related to the generator 
of displacements in the parameter 0, @| 



-2i[D UPt ], D t (B) 



dU t (B) r yt 



dB 



UJ(B). (3) 



The particular value of the Fisher information (and 
therefore SB) achieved in a given setting depends upon 
both the initial state of the atomic system po and the 
nature of the dynamical evolution Ut(B) [||. When ana- 
lyzing magnetometric limits, one typically takes the dy- 
namics to be generated by the Zeeman Hamiltonian 
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H(B) = -tvyBn-F, 



(4) 



FIG. 1: (color online) Schematic of a broadband atomic mag- 
netometer based on continuous observation of a polarized op- 
tical field double-passed through the atomic sample. 

where ft is the field unit vector B = Bri, 7 is the atomic 

gyromagnetic ratio and F a = YljLi /<* arc * ne c °l~ 
lective spin operators obtained from a symmetric sum 
over N identical spin-/ atoms (F = fN for a sample 
of N atoms each with total spin quantum number /). 
The time-evolution operator then satisfies ihd t Ut(B) = 
H(B)Ut(B) and optimizing the variance (A 2 D T (B)) over 
po for Eq. (H]) yields the shotnoise scaling 

SB T = l/ 7 rV2F (5) 

when po is separable, and the Heisenberg scaling 

SB T = 1/ 7 t2F (6) 

when entanglement is permitted between the different 
atoms. For N spin- 1/2 particles prepared into the initial 
cat-state (|ff • ■ ■ T) + l 44 '" 4))/v2 (in a basis set by ft), 
the uncertainty scaling under Eq. ((4]) is given by 1/jtN 
and often called the Heisenberg Limit. 

It was believed for some time that Eqs. ((5|) and (|6|) 
were fundamental: the 1/y/N scaling characteristic of 
the shotnoise uncertainty was an unavoidable byproduct 
of the coherent state projection noise (AF Z ) = \J F/2, 
only improved via entanglement; and 1/N scaling could 
only be achieved by creating a cat-state or maximally 
squeezed state, but never surpassed 0]. Recently, how- 
ever, it was shown that the 1/N "limit" can be overcome 
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[8j by extending the linear coupling that underlies E 
(U]) to allow for multi-body collective interactions 
Were one to engineer a probe Hamiltonian where B mul- 
tiplies fc-body probe operators, such as Fy, then the op- 
timal estimation uncertainty would scale more favorably, 
as ABk ~ 1/N k 8]. Unfortunately, the Zeeman Hamil- 
tonian is linear, meaning that l/N scaling is optimal for 
magnetometry based solely on Eq. (|4|). Improving mag- 
netometry beyond I/TV scaling therefore requires devel- 
oping some kind of nonlinear magnetic interaction. 

In this paper we provide evidence that one can improve 
upon the convencntional scalings in atomic magnetome- 
try, Eqs. (O and ©, by engineering effective dynamics 
that mimic a nonlinear coupling to the magnetic field 
using coherent positive feedback 2l|. Our approach in- 
volves double-pass ing an optical field through the atomic 
sample [HI U, 12, 3] at the same time that it is subject 
to Larmor precession, as depicted in Fig. (TT]). We support 
our claim of improved magnetometry by calculating the 
quantum Fisher information over a range of N sufficient 
to draw reasonable conclusions about the scaling. 

Effective Dynamics — To estimate a static magnetic field 
oriented along the y-axis B = By, the atomic sample is 
prepared into an ir-polarized spin coherent state. Like 
standard atomic magnetometers, a laser field with fre- 
quency uji , far-detuned from the nearest atomic resonance 
uj a , is used to observe Larmor precession by measuring 
the F z spin component via balanced polarimctry. On 
its first pass, the incoming linearly-polarized optical field 
acquires a Faraday rotation proportional to F z . Unlike 
standard Faraday spectroscopy, however, the probe laser 
passes through the atomic sample a second time [l(| prior 
to detection. Just before the second pass, a quarter wave- 
plate converts the first-pass Faraday rotation into eliptic- 
ity, which the atoms perceive as a fictitious magnetic field 
during the second pass. Since the second-pass field prop- 
agates parallel to the external magnetic field, the atoms 
see a total effective magnetic field that is the sum of B 
and the independent elipticity. As the atoms precess un- 
der B, the elipicity increases proportionately to F z and 
further amplifies F z (also the quantum uncertainty of the 
initial state), as illustrated by the simulation in Fig. 
which was generated using the model described below. 

To analyze the proposed magnetometer quantitatively, 
we assume that the interaction between the atoms and 
the far-detuned optical probe is well-described by a 
vector polarizability [3]. The collective atomic an- 
gular momentum F = (F^,F y ,F z ) couples to the 
two polarization modes of the traveling wave probe, 
which act as a Schwinger-Bose field with Stokes op- 
perators s = (sq, s x , s y , s z ). For convenience, we 
choose the polarization basis given by defining the time- 
domain Schwinger boson annihilation operator s t — 
\ J-^ 9{ u ) ®t uj&y ,.,s iut dui in terms of the annihilation 
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FIG. 2: (color online) Comparison of the double-pass mag- 
netometer depicted in Fig. ([Tj) to one using a conventional 
single-pass Faraday measurement. The simulated data was 
generated according to Eq. (I10|) with B = 0.1 (units of 7), 
F = 100, M — 1.0 and K = 1 X 10" 4 (units of 1/r). The 
single-pass comparison was obtained by setting K — 0. 



wave modes with frequency lu and form factor g(tu) ~ 
y^fku/AirceQ. The Stokes operators, s x = s + and 
s y = i(s — s^), are then reminiscent of quadrature op- 
erators and the interactions for each pass of the probe 
light through the sample can be expressed as = 
ih^F z (s\ - s t ) and = -hKF y {s\ + s t ) (Tl | . 



Developing a Markov approximation for the interac- 
tions, and h( 2 ', is a standard problem in the theory 
of open quantum svstems addressed by taking a weak- 
coupling limit [I5L fli| to obtain quantum Ito equations 



dUl 



(i) 



mF z {dS] - dS t ) - 2 mF z dt f u t 



du[ 2) = S.iVkFy(dSl +dS t )-^kF^dt^ul 2) 



(7) 



(8) 



for the first- and second-pass unitary propagators [11 
12 1 . In these expressions, dS\ and dS t are delta- 
correlated noise operators derived from the quantum 
Brownian motion St = J s u du. The noise terms sat- 
isfy the quantum Ito rules: dStdSj = dt and dSjdSt = 
dSf = {dS\f = 0, and can be viewed heuristically as 
a consequence of vacuum fluctuations in the probe field. 
The rates m = 2ir\fig(uji)\ 2 and k — 2ir\Kg(oji)\ 2 can be 
computed using the methods found for example in Rcfs. 
[l6| . Assuming that the time-scale for the the light to 
propagate (twice) through the same is fast compared to 
the magnetic dynamics, the separate evolutions, Eqs. ([7]) 
and ([5]), can be combined into a single Markov limit, 
yielding the propagator for the double-passed dynamics 



dU t = 



M 



K 



iVMKF y F z --H- —Fl - -F 2 ) dt (9) 



MF z (dS J t -dSt) 



KFy{dSt+dS t ) U u 



for the x and y-polarzed plane- 



where M = m/2, K = k/2 and H = H(B) is the Zeeman 
Hamiltonian from Eq. ([JJ [ll| . 
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The measurement used to estimate B is obtained by 
detecting the optical helicity of the probe laser, as this 
observable caries information about the z-component 
of the atomic spin. A continuous measurement of the 
helicity is described by the time-evolved Stokes oper- 
ator Z t — U}{S\ + St)Ut, with the dynamics dZ t — 
2VMU}P z U t dt + {dS} + dS t ). It can be shown that 
the observations constitute a classical stochastic pro- 
Zt 16, The techniques of quantum filter- 

17, Il8| can then be used to obtain the 
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cess 

ing theory 

best least-squares estimate of the expectation value of 
any atomic operator X conditioned on the measurement 
TTt(X) = ]&[X\Zt] = tr[Xpt] where the conditional den- 
sity operator satisfies 



dp t = i 7 B[F y ,p t }dt + VKM[F y ,{F z ,pt}]dt 

+MV[F z ]p t dt + KT>[F y }p t dt (10) 
+ (VMM[F z ]pt + iVK[F y ,pt])dW t . 

Here, the innovations dW± = dZ t — 2y/Mtr[F z pt]dt are a 
Wiener process, E[tWt] = 0, dW? — dt, and the super- 
operators are defined as V\X]p t = Xp t X^ — (X'Xpt + 
p t X^X)/2 andM[X]p t = Xpt + p t X-tT[{X + X^pt]pt. 
A derivation of Eq. (fT0|) can be found in Ref . [ll| . 



Quantum Fisher Information Results — The quantum fil- 
ter Eq. (|10p can be used to simulate individual measure- 
ment realizations of the double-pass magnetometer [c.f., 
Fig. ((21)] • Furthermore, one can evaluate the quantum 
Fisher information (at least numerically) via a finite- 
difference approximation to the logarithmic derivative 



dp t (B) p t {B + dB) - p t (B - dB) 



dB 



dB 



(11) 



by evolving pt(B), p t (B + dB) and Pt(B — dB) under the 
same noise realization for dB <C 1 [Hi . On an individual 
trajectory basis, the Fisher information calculated using 
the conditional density operator T t \Z t = tr[£fpt(-B)] is 
also conditioned on the particular measurement realiza- 
tion that generated pt and must be averaged over many 
realizations to obtain the unconditional quantum Fisher 
information T t = TE[T t \Zt\. The lower bound SB T can 
then be obtained from Eq. (pj with statistical errorbars 
given by a(SB r ) = ^' 2 a[l T \Z t }/2. 

We calculated 1t(B) over a range of spin quantum 
quantum numbers F — Nf spanning more than an order 
of magnitude to determine a lower bound on the mag- 
netic field estimation uncertainty using Eq. ([1]). Our re- 
sults indicate that the Fisher information depends heav- 
ily upon the choice of the coupling strengths M and K, 
which is not surprising since the measurement strength 
M determines how much spin-squeezing is generated and 
K determines the strength of the effective nonlinearity. 
Like any measurement procedure that involves amplifica- 
tion, both the signal and noise are affected, and optimal 
performance requires choosing the correct gain. 
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FIG. 3: (color online) Comparison of the estimation uncer- 
tainty AB as a function of the total atomic angular mo- 
mentum (proportional to iV) for double-pass and single-pass 
atomic magnetometers determined by calculating the quan- 
tum Fisher Information with M = 1 (in units of 1/r) and 
K = 1 x 10~ 4 chosen to be optimal for F = 140/i . 



If one choses M = 1/r, to obtain an optimal spin- 
squeezed state at the final time t = r [5|, then it is 
straightforward to optimize over the nonlinearity K, as 
illustrated in the inset of Fig. © for F = \Q0h. We 
found that the optimal value K * depends upon the num- 
ber of atoms, and that the Fisher information saturates 
and then decreases if the number of atoms exceeds the 
value of N — F sat /f used to compute K*(F). Figure 
[3]) shows the behavior of SB T as a function of F up to 
the saturation point F < F sa t ~ 150. The largest value 
of F prior to saturation yields a 8B T that is slightly be- 
low the bound 1/tjF 3 ^ 2 that would be obtained for a 
two-body coupling Hamiltonian and an initially separa- 
ble state po [19j . Despite this saturation of the quantum 
Fisher information for F > F sat at a given choice of K, 
one can choose the value of K* such that saturation oc- 
curs only for i^ sat > Flnax over any specified finite range 
F < -Fmax- An improvement beyond 1/N scaling can be 
achieved over any physically realistic number of particles. 

A second approach to avoiding saturation of the Fisher 
information for large F is to scale the parameters M and 
if as a decreasing function of F. For practical considera- 
tions, it is also desirable to set M = K as these parame- 
ters are determined by the atom-field coupling strengths 
on the first and second pass interactions, thus quanti- 
ties such as the laser intensity and detuning not easily 
changed between the two passes. We have found that 
scaling M and K according to the functional form 



M =K = c/rF a , 



(12) 



where c and a are constants, leads to a power-law scaling 
for the uncertainty bound SB T ~ 1/N k . The inset plot 
in Fig. (01 shows the slope of a linear fit of log 10 5B T to 
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FIG. 4: (color online) Evidence that the field estimation un- 
certainty AB can be made to scale as a power law SB T ~ 
1/N by decreasing the parameters M and K as a function 
of the total angular momentum F according to Eq. (|12|) with 
ttfx3/4. The power-law fit (solid line) has a slope of —0.97. 

log 10 F (i.e., a slope of k = — 1 corresponds to the Heisen- 
berg uncertainty scaling) as a function of a (with c chosen 
so as to avoid the saturation behavior described above). 
As demonstrated by the data points in Fig. (TJ|, it is pos- 
sible to achieve 1/N scaling (to within a small prefactor 
offset) with a = 0.77 and c = 0.589. The distribution 
of conditional uncertainties SB T \ Z t for the statistical en- 
semble of measurement realizations [dots in Fig. ([JJ] is 
depicted for the different values of F. The mean and un- 
certainty of this distribution are denoted by the circles 
and errorbars, and a fit to this data gives 8B T ~ ^-0.97^ 
Extrapolating the values of M and K that would be 
required to achieve a sensitivity SB T corresponding to 
the Heisenberg uncertainty scaling Eq. ([6]) with TV ~ 10 8 
(small for typical experiments) implies that one would 
only require Mr ~ 10~ 7 . Since Mr = 1 corresponds 
to a maximally spin-squeezed state at time t = r, these 
results suggest that 1/N scaling can be achieved by an 
extremely weak measurement, with no appreciable gener- 
ation of conditional spin-squeezing. These results suggest 
that it may be much easier than previously believed to 
to outperform the shotnoise scaling since doing so should 
not require generating an entangled state of the atoms. 

Conclusion — Our results highlight that there are two 
complementary ways to improve metrological sensitivity: 
(1) reducing the quantum noise of the probe; and (2) 
enhancing the accumulated phase acquired by the probe 
due to its interaction with the metrological quantity. Our 
procedure does both — probe uncertainty is reduced by 
conditional spin squeezing, and phase accumulation is en- 
hanced by coherent optical feedback. It does, however, 
remain an open question to reconcile precisely how the 



effective dynamics observed here can be understood using 
the Hamiltonian framework of Ref. |8j. This outstand- 
ing problem is extremely important since effective multi- 
body interactions, like those developed here, are likely 
to be essential to optimizing sensitivity in precision mea- 
surements. This work was supported by the NSF (PHY- 
06399 94) and the AFO SR (FA9550-06-0 1-0178). Please 
visit http://qmc.phys.unm.edu/ to download the simula- 
tion code used to generate our results as well as all data 
files used to generate the figures in this paper. 
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